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The s process in massive stars at low metallicity. Effect of 
primary 14 N from fast rotating stars. 
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ABSTRACT 

The goal of this paper is to analyze the impact of a primary neutron source on 
the s— process nucleosynthesis in massive stars at halo metallicity. Recent stellar 
models including rotation at very low metallicity predict a strong production 
of primary 14 N. Part of the nitrogen produced in the H— burning shell diffuses 
by rotational mixing into the He core where it is converted to 22 Ne providing 
additional neutrons for the s process. We present nucleosynthesis calculations 
for a 25 M star at [Fe/H] = —3, —4, where in the convective core He— burning 
about 0.8 % in mass is made of primary 22 Ne. The usual weak s— process shape is 
changed by the additional neutron source with a peak between Sr and Ba, where 
the s— process yields increase by orders of magnitude with respect to the yields 
obtained without rotation. Iron seeds are fully consumed and the maximum 
production of Sr, Y and Zr is reached. On the other hand, the s— process efficiency 
beyond Sr and the ratio Sr/Ba are strongly affected by the amount of 22 Ne and by 
nuclear uncertainties, first of all by the 22 Ne(a,n) 25 Mg reaction. Finally, assuming 
that 22 Ne is primary in the considered metallicity range, the s— process efficiency 
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decreases with metallicity due to the effect of the major neutron poisons 25 Mg 
and 22 Ne. This work represents a first step towards the study of primary neutron 
source effect in fast rotating massive stars, and its implications are discussed in 
the light of spectroscopic observations of heavy elements at halo metallicity. 

Subject headings: stars: abundances - - chemically peculiar - - early-type - 
rotation 



Introduction 



Very metal— poo r halo stars have form ed from matter that was only enriched by mas- 
sive stars ejecta (e.g. iTravaglio et al.ll2004 , hereafter Tr04). Their surface compositions thus 
reflect the nucleosynthesis occurring in the first generations of massive stars. A puzzling 
feature observed in a si gnificant number o f stars is the high Sr enrichment correlated with 



low Ba abundance (e.g. Truran et al. 2002 ; Aoki et aDbood ). with the observed Sr/Ba ratio 



spreading over more than two orders of magnitude. They are not explained by the standard 
r process (rapid neutron capture process. iBurbidge et al.lll957l ). In particular, the r— process 
is responsible for t he production of a bout 20 % of solar Ba and 12 % of solar Sr, causing a 
[Sr/Ba] r = -0.22 JSneden et alil2008h . On the other hand, stars like CS 22949-037 ([Fe/H] 
= —4) and the most metal— poor star known HE 1 327—2326 ([Fe/H] = —5.45) are Sr— Y 



rich, Ba poor and with high [Sr/Ba] (for details see iDepagne et al.l 120021 ; lAoki et ajj [2006 



respectively). Another process must be responsible for the Sr observations, w hich is related 
to a similar observed Y and Zr overproduction with respect to Ba (e.g., Tr04. iMontes et al. 



20071 ). Asymptotic giant branch (AGB) stars do not have time to contribute to the Galac- 
tic chemical evolution (GCE) at very low metallicity. The standard weak s process (slow 
neutron capture process) in massive stars can neither be the solution for the Sr— Y— Zr over- 
production. Indeed, the secondary nature @ of the main neutron source 22 Ne progressively 
reduces the s— process efficiency with decreasing metallicity. In particular, for [Fe/H] < —2 
primary neutron poisons capture most of the neutrons causing an even faster decrease of the 
s— process yields with decreasing metallicity (IRaiteri et al.lll992l ). Tr04 propose that an un- 
known primary neutron capture process in massive stars is responsible for the production of 
Sr— Y— Zr, and provide an evaluation of the amount of Sr, Y and Zr that should be produced 
in order to explain observations at very low metallicity: about 10 % of the solar Sr and 20 % 
of the solar Y and Zr. They called this process light element primary process (LEPP). The 



8 The production of secondary (primary) isotopes does (does not) depend on the initial metallicity of the 
star. 
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unknown component could also be produced by explosive nucleosynthesis in massive stars. 
If in core collapse supernovae the a— rich freezout is not complete, the final explosive nucle- 
osynthesis in high— entrop y neutrino winds from a fo rming neutron st ar may be significant 
for t he Sr peak elements ( Wooslev fc Hoffman I fl992l ) . More recently, Farouqi et al. ( 2008 ) 
and iQian fc Wasserburg I (120071 ) generally confirm the feasibility of this nucleosynthesis sce- 
nario in high— entropy neutrino winds. 

Recent results from massive stars calculations at low metallicity including rotation could 
have important implications in this discussion. In massiv e stars at very low m etall icity, ro 



tation boosts the formation of primary 14 N. According to lMeynet et al.l (120061 ) and iHirschi 
( 120071 ) . due to rotational mixing primary 12 C produced in the convective He core is brought 
in contact with the bottom of the convective H shell, where it is converted into 14 N via the 
CN cycle. P art of the 14 N is mixe d in the envelope, and par t is mixed down or en gulfed in the 
He core (e.g. iMeynet et al.ll2006l ). A more recent model by lHirschi et al.l (120081 ) for a 20 M Q 
star predicts that H shell ashes are directly processed by the convective core He— burning, 
causing an even higher amount of primary 14 N in the He core than in the previous mod- 
els. During He— burning, nitrogen is rapidly converted to 22 Ne by the a— capture sequence 
14 N(a,7) 18 F(/3+i/) 18 0(a,7) 22 Ne. Close to the He exhau stion, 22 Ne(a,n) 25 Mg is activated pro- 
viding neutrons for the s process. IHirschi et al.l (120081 ) quantified that for a 20 M Q star with 
[Fe/H] = —4.3 the primary 22 Ne in the He core before the s— process activation is about 0.8 
% in mass, about 200 times m ore than in t he cas e wit h no rotation. 

Spectroscopic observations by ISpite et al.l (120051 ) and llsraelian et al.l (120041 ) indicate a pri- 
mary production of nitrogen over a large metallicity range for [Fe/H] > —4.1, and also at 
lower metallicities strong N enr ichment has been o bserved in stars c orrela ted with Sr en- 
richment (e.g. HE 1327-2326. iFrebel et all boosh . Ichiappini eiTall (bood ) show that the 
primary nature of nitrogen can be explained at low metallicity by using the yields of fast 
rotating massive stars. Since the production of primary nitrogen is directly related to the 
22 Ne available in the He core, we may also expect that the additional 22 Ne is produced in a 
primary way over a large metallicity range. 

In the present work, our aim is to investigate possible consequences of massive star models 
including rotation on the synthesis of s— process elements. In Sj2] we present s— process nu- 
cleosynthesis calculations in the He core and in the following convective C— burning shell of 
a 25 M at [Fe/H] = —3, —4. Final discussion and conclusions are given in §3] 



2. Calculations and results 



We present s— process calculations for a 25 M star at [Fe/H]= —3, —4. A post— processing 
code has been used to follow the s process in the convective He core phase and in the convec- 
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tive C shell phase ( IRaiteri et al.lll99ll ). In particular, the C shell burns over the ashes of the 
previous He core. In this phase the main neutron source for neutron capture is the 22 Ne left 
after He— burning. During the supernova explosion, most of the convective C shell region is 
ejected unchanged carrying the s— process signature of the pre— explosive phase. For a 25 
M , indeed, the material up to 3—3.5 M is further processed and the s yields are destroyed. 
Between 3—3.5 M and 6—6.5 M the O rich arid s— process rich material is ejected almost 
unchanged by the explosion (e.g. IWoosley et al.ll2002l ). 

The nuclear network has been updated, as described in iPignatari et al.l (120061) . In par- 
ticul ar, the 22 Ne(a,n ) 25 Mg and the 12 C(a,7) 16 rates are from Jaeger et al.l ( 2001 ) and 
from ICaughlan et al.l (119851 ). respectively. I n the initial compositio n, a— enhancements for 
the light isotopes are i nclud ed ac cording to iMishenina et al.l ( 120001 ) for the oxygen and to 



Goswami &: Prantzos I (120001 ) and iFrancois et al.l (120041 ) for the other elements. Heavy iso- 



tope abundances beyond iron are solar— scaled with the initial metallicity. 
In Fig. [1] Top Panel the isotopic distributions between 57 Fe and 138 Ba are presented at the 
end of the C— burning shell of a 25 M star at [Fe/H] = —4 for three cases: i) the standard 
s— process distribution, ii, Hi) s— process distributions where the additional 22 Ne is 0.5 % 
and 1.0 % in mass (about 100 and 200 times more than in the standard case), respectively, 
in agreement with iHirschi et al.l (120081 ) . For each of these cases, in Table Q] the s— process 
parameters have been reported with the final C shell abundances Y; for a sample of elements. 
The ejected C shell mass of element i, em;, is given by em; = Dm x A; x Yj, where Ai is 
the average mass number of element i and it is given in Table (H and Dm is the C shell mass 
ejected unchanged by the explosion in M . For comparison, the amount of 56 Fe produced 
by a 25 M star is 0.1—0.2 M , it i s primary and it strongly depends on the explosion 
parameters (e.g.. iRauscher et al.ll2002l ). 

In case i the isotopic distribution is peaked between Fe and Sr. The 22 Ne abundance depends 
on the initial CNO abundances used for a given metallicity and on the a— enhancements used 
for light elements. As already mentioned in SJTJ for [Fe/H] < —2 the main neutron poisons are 
primary in the He core, and the probability that a free neutron is captured by iron seeds is so 
low that the s— process efficiency de creases more rapidly than the iron seeds with decreasing 



metallicity (e.g. IRaiteri et al.lll992l ). Notice that in the C shell the s— process efficiency is 
strongly reduced by primary neutron poisons like 20 Ne and 23 Na (produced by the 12 C + 12 C 
reaction). At [Fe/H] = —4, in cases ii and Hi iron seeds are more consumed with respect to 
standard s— process calculations. The s— process efficiency between Sr and Ba is increased 
by orders of magnitude, producing a flat isotopic pattern in this mass region. In particular, 
case ii is more efficient in producing isotopes in the Pd region and less efficient in producing 
heavier elements with respect to case Hi. On the other hand, case ii and case Hi provide 
similar results for Sr, Y and Zr (within a factor 1.3), at the bottle— neck of the neutron 
capture flow. This implies that the Sr/Ba ratio decreases with increasing 22 Ne: according 
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to Table [H in case ii the Sr/Ba = 9.1 and in case Hi Sr/Ba = 3.8. The addition of more 
primary 22 Ne does not imply a higher production of Sr, Y and Zr. This happens for two main 
reasons. The first reason is related to the light isotope abundances. Indeed, at [Fe/H] = —4 
most of the neutrons produced by 22 Ne are captured by the strongest neutron poisons, 25 Mg 
and 22 Ne ( 16 is the strongest neutron absorber, but the neutrons are efficiently recycled 
via the 17 O(a,n) 20 Ne reaction). The second reason is related to the s— process seeds. In Fig. 
[I] Top Panel Fe seeds have been destroyed by the s— process, and the resulting abundances 
for elements between Fe and Sr are smaller than the abundances of elements beyond Sr. In 
case ii of Fig. [1] Top Panel the resulting s— process production yields for elements up to 
Rb (included) is less than 1 % than for elements beyond it. This implies that the neutron 
capture nucleosynthesis flow producing Sr is exhausted. 



In Fig. [T] Bottom Panel we present similar results at [Fe/H] = —3. The correspondent 
s— process parameters and element abundances h ave b een reported in Table [H The new 
generation of models presented by iHirschi et al.l (120081 ) does not yet include models with 
[Fe/H] ~ —3. However, since the nitrogen is observed to be primary over a large range 
of metallicity (ISpite et al.l 120051 ) . we discuss now the possibility that the additional 22 Ne is 
produced with similar efficiency in the He core at [Fe/H] = —3. In Fig. [1] Bottom Panel the 
isotopic distribution wit hout additional 22 Ne ( case i) is similar to t he isotopic distribution 
at solar metallicity (e.g. iRauscher et al.l 120021 : iPignatari et al.ll2006l ). In particular, a peak 
of production is obtained in a large mass region between Cu and Sr and with a rapid fall 
beyond Y, where the standard weak s process is not efficient. On the other hand, the 
production factors in case i are about three orders of magnitude lower than in the case at 
solar metallicity, confirming the secondary nature of the weak s— process in massive stars. 
In case ii and case Hi the primary 22 Ne is 0.2 % and 1 % in mass (about 10 and 50 times 
more than in the standard case with no rotation), respectively. In this case we consider 
a larger range of 22 Ne with respect to [Fe/H] = —4, to test the sensitivity of s— process 
calculations to different amounts of the neutron source. Like at [Fe/H] = —4, the mass 
region Fe— Zn is again substantially depleted in cases ii and Hi with respect to case i. On 
the other hand, quite similar results are observed between Zn and Kr. Finally, abundances 
beyond Kr are strongly enhanced in the calculations including primary 22 Ne. At [Fe/H] = 
—3 the iron seeds are a factor ten higher than at [Fe/H] = —4, and the s— process efficiency 
increases by about a factor ten in cases ii and Hi with respect to similar cases in Fig. [1] Top 
Panel. As recalled above, since the neutron source 22 Ne and the main neutron poisons 25 Mg 
and 22 Ne are primary, the s— process is secondary scaling with the iron seed abundance. 
Without considering smaller local differences, the isotopic patterns at [Fe/H] = —3,-4 are 
quite similar. Notice that changing the primary 22 Ne amount by a factor five affects the Sr, 
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Y and Zr production less than a factor two. Concerning heavier isotopes, more primary 22 Ne 
implies a more efficient production in the Ba region, more 22 Ne left in the He core ashes and 
as a consequence higher neutron densities in the C shell (e.g., 96 Zr overproduction in case 
iii with respect to case ii). In particular, in case ii Sr/Ba = 107.6, which is strongly larger 
than in case iii with Sr/Ba = 4.6. 

Considering metallicities higher than [Fe/H] ~ —3 and closer to solar metallicity, rotational 
velocity and the related mixing efficiency in stars should decrease, assumin g that stars of 



different metallicity begin their evolution with similar angular momentum (e.g. lMeynet et al. 



20061 ; iHirschi 1 120071 ). Therefore, the production of primary 14 N and of the additional 22 Ne 



in the He core becomes negligible at solar metallicity. This implies that the new s— process 
component discussed in this paper is activated at metallicities typical of the Galactic halo, 
but not for metallicities of the Galactic disk. 

An interesting point to discuss is how much the nuclear rates affect our results. It is well 
known that the 22 Ne(a,n) 25 Mg reaction rate and the relative (a,n)/(a,7) ratio affect the 
s— process nucleosynthesis. In Fig. [2] we present case ii i of Fig. [1] Bottom Panel where 
the 22 Ne(ci!,n) 25 Mg recommended rate (j Jaeger et al.ll200ll ) has been multiplied and divided 
by a factor two. Elements lighter than Sr are more produced with a lower 22 Ne(a,n) 25 Mg 
rate, since less neutrons are available to feed heavier elements through following neutron 
captures. The opposite effect is observed with higher 22 Ne(o;,n) 25 Mg. Higher (a,n) rate 
and higher (o;,n)/(a,7) ratio do not strongly affect our calculations between Sr and Ba. 
The reason is that in this case more neutrons are produced but also more 25 Mg, which is 
the strongest neutron poison in these conditions. Furthermore, iron seed are destroyed and 
cannot feed the nucleosynthesis of Sr, Y and Zr. For instance, in case iii of Fig. [T] Bottom 
Panel the amount of material between Fe and Rb is less than 1 % of the material beyond 
Sr. On the other hand, the s— process distribution beyond Sr shows a strong propagation 
effect due to a lower (a,n) rate. In this case, less neutrons and 25 Mg are produced but more 
22 Ne is left, which is also a neutron poison, causing a lower s— process efficiency. Assuming a 
factor two of uncertainty for the 22 Ne(a,n) 25 Mg, the Sr/Ba ratio changes by a factor 221 (0.5 
and 110.4 for the (ct,n)*2 case and for the (a,n)/2 case, respectively)! Out of the range of 
Fig. [1] Bottom Panel, the distribution beyond Ba is also significantly affected. For example, 
the ratio Sr/Pb ranges between 1.9 and 8059.7 for the (a,n)*2 case and the (a,n)/2 case, 
respectively. This large effect is due to the amount of 22 Ne consumed in the He core: Sr, Y 
and Zr are no more produced, but heavier elements in Ba and Pb peaks are fed at the expense 
of lighter s— process isotopes. Related to this last point is also the nuclear uncertainty of 
12 C(a,7) 16 0, which is in competition with the 22 Ne(a,n) 25 Mg capturing helium before its 
exhaustion in the He core. For instance, a lower 12 C(a,7) 16 rate increases the amount of 
22 Ne burnt in the He core, and it implies a more efficient s nucleosynthesis at both the Ba 
and the Pb peaks. Other nuclear uncertainties (e.g. the 16 0(n,7) 17 cross section) show 



- 7- 



smaller but non— negligible effects on the calculations. 



3. Discussion and final remarks 

We presented a first study of the effect of rotational mixing in massive stars on the s 
process at [Fe/H] = —3,-4. A grid of models for a large spread of masses and metallicities 
is required for a more detailed analysis, including an evaluation of the impact of stellar 
model physics uncertainties. Notice also that other processes not considered in the present 
models may affect the mixing efficie ncy and the nucleos ynthesis in massive stars, e.g., binary 



interaction and magnetic field (e.g. lLanger et al.ll2008l ) 



In the calculations presented at [Fe/H] = —3,-4, the additional 22 Ne due to rotation is 
assumed to be primary in the He core and in the following C shell, and it causes a peak of 
production in the region between Sr and Ba. In this mass region the s— process efficiency 
is orders of magnitude higher than in the case with no rotation. The main neutron poisons 
are 25 Mg and 22 Ne. Iron seeds are strongly depleted feeding heavier isotopes, and once 
they are destroyed the production flux to Sr— Y— Zr is substantially exhausted. In these 
conditions, the Sr— Y— Zr s nucleosynthesis is not more efficient increasing the 22 Ne burnt 
in the He core, but Sr— Y— Zr elements can feed the production of heavier elements. The 
s— process efficiency at the Ba peak (and at the Pb peak) is affected by the amount of 
primary 22 Ne, by the fraction of 22 Ne burnt in the He core and by nuclear uncertainties, in 
particular by the 22 Ne(a,n) 25 Mg and the 12 C(a,7) 16 0. Indeed, the Sr/Ba ratio changes by 
two orders of magnitude considering a factor of two of uncertainty for the 22 Ne(a,n) 25 Mg. 
The s process in fast rotating massive stars could provide a new scenario to explain the Sr 



enric hment coupled with high Sr/Ba observed in many stars (ITruran et al.ll2002l ; lAoki et al. 



20051 ) at very low metallicity, including also the spread in the observed Sr/Ba ratio. On 
the other hand, the s— process component activated by the additional 22 Ne due to rotation 
is a secondary— like process for the Sr peak elements in the range of metallicity considered, 
whereas the Sr enrichment is still observed in stars like HE 1327—2326 ([Fe/H] = —5.45). 
This last point could be a strong constraint to rule out this scenario. Unfortunately, elements 
between Zr and Ba have been observed in few metal poor stars to constrain the signature 
of the proposed s— process mechanism. A noticeable exception is the star HD 122563 (e.g., 



Montes et al.l 120071 . and references therein), where one— two lines have been detected and 
measured for elements in this mass region. More observations are required to draw firm 
conclusions. New models at higher metallicity are necessary to study the chemical evolution 
of the Galactic disk. Models presented here could have important consequences for the 
chemical evolution of elements between Sr and Ba in the Galactic halo. 
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Table 1: For the cases in Fig. [I] the neutron exposure r, the peak central neutron density in 
the He core and the peak neutron density in the C shell are reported. The initial 22 Ne mass 
fraction available in the He core and the final C shell number abundances Y; for a sample 
of elements are reported, and in the first column the average mass number of element i is 
given (Aj). 



M = 25 M 


[Fe/H] = -3 


[Fe/H] = -3 


[Fe/H] = -3 


[Fe/H] = -4 


[Fe/H] = -4 


[Fe/H] = -4 


Y( 22 Ne) ini 


2.12 x 10" 4 


2.0 x 10 -3 


1.0 x 10 -2 


5.21 x 10~ 5 


5.0 x 10 -3 


1.0 x 10 -2 


r (mb -1 ) 


0.255 


0.659 


0.923 


0.118 


0.841 


0.942 


Heco(n n ) max (cm -3 ) 


2.808 x 10 7 


6.073 x 10 7 


7.912 x 10 7 


1.475 x 10 7 


7.317 x 10 7 


7.976 x 10 7 


Csh(n n ) max (cm -3 ) 


7.278 x 10 11 


1.518 x 10 12 


3.691 x 10 12 


7.021 x 10 11 


2.969 x 10 12 


4.362 x 10 12 


Sr (87.7) 


9.304 x 10 -11 


7.440 x 10 -9 


5.595 x 10~ 9 


5.823 x 10 -13 


6.589 x 10 -10 


5.294 x 10 -10 


Y (89.0) 


6.964 x 10" 12 


1.746 x 10 -9 


1.992 x 10 -9 


6.450 x 10 -14 


2.130 x 10 -10 


1.935 x 10 -10 


Zr (91.3) 


5.754 x 10 -12 


2.806 x 10 -9 


5.094 x 10 -9 


6.623 x 10 -14 


4.745 x 10 -10 


5.115 x 10 -10 


Mo (96.0) 


3.751 x 10 -13 


3.656 x 10 -10 


6.885 x 10 -10 


6.663 x 10 -15 


1.055 x 10 -10 


1.526 x 10 -10 


Ba (137.4) 


5.222 x 10" 13 


6.917 x 10 -11 


1.222 x 10" 9 


3.967 x lO -14 


7.270 x 10 -11 


1.378 x 10" 10 
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Fig. 1. — TopPanel: s— process distributions between Fe and 138 Ba normalized to solar for 
the 25 M Q and [Fe/H] = —4 at the end of the convective C— burning shell. The horizontal 
lines correspond to the 16 overabundance in the C shell (thick line), multiplied and divided 
by two (thin lines). Isotopes of the same element are connected by a line. The cases presented 
are the following: i) standard (black triangles), X( 22 Ne) ini = 5.21 x 10~ 5 ; ii) X( 22 Ne) ini = 
5.0 x 1CT 3 (open squares, full squares for the s— only isotopes); Hi) X( 22 Ne) ini = 1.0 x 10~ 2 
(open circles, full circles for the s— only isotopes). BottomPanel: the same as Top Panel, 
at [Fe/H] = -3: i) standard (black triangles), X( 22 Ne) ini = 2.12 x 10~ 4 ; ii) X( 22 Ne) ini = 
2.0 x 10~ 3 (small triangles); Hi) X( 22 Ne) ini = 1.0 x 10~ 2 (small circles). 
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Fig. 2. — The same as case Hi in Fig. [T] Bottom Panel, but using different 22 Ne(a,n) 25 Mg 
rates: recommended (circles), recommended * 2 (squares) and recommended / 2 (triangles). 



